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Abstract 

Gravity modes are the best probes to infer the properties of the solar radiative 
zone that represents 98% of the Sun's total mass. It is usually assumed that high- 
frequency g modes give information about the structure of the solar interior whereas 
low-frequency g modes are more sensitive to the solar dynamics (the internal ro- 
tation). In this work, we develop a new methodology, based on the analysis of the 
almost constant separation of the dipole gravity modes, to introduce new constraints 
on the solar models. To validate this analysis procedure, several solar models - in- 
cluding different physical processes and either old or new chemical abundances (from, 
respectively, Grevesse and Noels (Origin and evolution of the elements. Cambridge, 
England, 199, 15, 1993) and Asplund, Grevesse, and Sauval (Astron. Soc. Pac. CS, 
San Francisco, 336, 25, 2005)) - have been compared to another model used as a 
reference. The analysis clearly shows that this methodology has enough sensitivity to 
distinguish between some of the models, in particular, between those with different 
compositions. The comparison of the models with the g-mode asymptotic signature 
detected in GOLF data favors the ones with old abundances. Therefore, the physics 
of the core - through the analysis of the g-mode properties - is in agreement with 
the results obtained in the previous studies based on the acoustic modes, which are 
mostly sensitive to more external layers of the Sun. 
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1. Introduction 

The combination of detailed modeling and helioseismic observations has enabled us 
to accurately describe the solar interior from the surface to its center where the 
nuclear reactions take place. The precise characterization of the resonant acoustic 
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(pressure or p) modes allowed constraining the theories describing the solar interior. 
The profiles of sound speed, density (e.g. Basu et at, 1997; Couvidat, Turck-Chieze, 
and Kosovichev, 2003), and differential rotation (Thompson et al, 1996; Howe et al., 
2000; Chaplin et al, 2001; Couvidat et al., 2003), as well as the position of the base 
of the convective zone (e.g. Christensen-Dalsgaard, Gough, and Thompson, 1991), 
and the prediction of the neutrino fluxes (e.g. Turck-Chieze et al, 2001) are some 
examples of the constraints provided by the p-mode analyses. 

Unfortunately, only a few p modes propagate inside the inner core - below 0.25 
solar radius (Rq) - containing about 50% of the solar mass. To progress in the 
knowledge of these layers, another type of mode is needed: the gravity (g) modes. 
These modes have not yet been unambiguously detected individually in the Sun 
(for example see Appourchaux et al, 2000; Gabriel et al, 2002; Elsworth et al, 

2006) because they become evanescent in the convection zone reaching the surface of 
the Sun with very small amplitudes (Andersen, 1996; Kumar, Quataert, and Bahcall, 
1996). Some candidates have been detected (Turck-Chieze et at, 2004; Mathur et al, 

2007) analyzing helioseismic data from the Global Oscillations at Low Frequencies 
(GOLF (Gabriel et al, 1995)) instrument on board the ESA/NASA Solar and Helio- 
spheric Observatory (SOHO (Domingo, Fleck, and Poland, 1995)) spacecraft. These 
candidates have provided some scenarios about the possible dynamics of the solar 
core but it has been impossible to correctly label these peaks in terms of their orders 
and degrees (n, £, m) and to give the central frequencies of the candidates. Recently 
Garcia et al. (2007) were able to uncover the signature of the asymptotic properties of 
the dipole (£ = 1) g modes measuring their constant spacing APi. Indeed, for a given 
degree I, the difference APi of the periods of g modes with consecutive radial order n 
is almost constant when n»< (Tassoul, 1980; Provost and Berthomieu, 1986). The 
comparison of the observations with some solar models - including different rotation 
rates of the core - showed the great sensitivity of this analysis to the dynamics of 
the solar core and it implied a rotation rate in the core faster than in the rest of the 
radiative region. 

During the last few decades, many new physical processes have been studied and 
included in the solar models while others have been refined in order to better match 
the helioseismic observations. Today, these models include better descriptions of the 
microscopic physics such as opacities, equation of state, nuclear reaction rates (e.g. 
Morel, Provost, and Berthomieu, 1997 and references therein), microscopic diffusion 
(Michaud and Proffitt, 1993) as well as other processes, for example turbulence 
mixing in the tachocline (Brun, Turck-Chieze, and Zahn, 1999). 

New studies of the solar spectrum using the most recent advances in 3D mod- 
eling of the solar atmosphere and taking into account Non Local Thermodynamic 
Equilibrium (NLTE) effects have provided new estimations of the solar photospheric 
composition (Asplund, Grevesse, and Sauval, 2005). In comparison with the older 
estimations (e.g. Grevesse and Noels, 1993) based on ID models and the LTE approx- 
imation, this new measurement leads to a smaller metallicity and to a substantially 
modification of the abundances of C, N, and O, which are crucial in solar modeling. 
The standard models computed with these new abundances yield a larger discrepancy 
with the helioseismic observations. They disagree with sound-speed and density 
profiles inferred from acoustic modes (for example see Turck-Chieze et al, 2004; 
Bahcall, Serenelli, and Basu, 2005; Guzik, Watson, and Cox, 2005) and we need, for 
example, to modify some opacities to reduce the discrepancy. It has also been shown 
that the models computed with these new abundances modify the frequencies of the 
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resonant low-degree p and g modes (Basu et al, 2007; Chaplin et al., 2007; Zaatri 
et al, 2007; Mathur et al, 2007). 

In this paper we use a set of solar models (briefly described in Section 2) - 
including different physical processes combined with the old and new abundances of 
the Sun and with different core rotation rates - that we utilize as a tool to study 
whether the methodology used to detect the dipole gravity modes in the Sun by 
Garcia et al. (2007) (and summarized in Section 3) is able to distinguish between 
these models. To do so, we compare in Section 4 all of them to our reference model: 
the Saclay seismic model (Turck-Chieze et al, 2001). Finally, once the method has 
been demonstrated, we compare, in Section 5, the solar models with the GOLF data. 



2. Description of the Solar Models 

Several solar models have been computed with the Code d'Evolution Stellaire Adap- 
tatif et Modulaire (CESAM) (Morel, 1997), which differ from each other by the 
physical processes that they include. They have been developed to study the influence 
of different physical processes and the impact of the new abundances on the (/-mode 
frequencies. These models are fully explained in Mathur et al. (2007). Here we recall 
briefly their main characteristics. 

First let us describe the solar model that has been taken as a reference: the 
so-called seismic model (corresponding to the model seismic^ of Couvidat, Turck- 
Chieze, and Kosovichev, 2003). It is important to note that this model is produced 
by an evolutionnary code and is not directly deduced from a seismic inversion. It 
is "seismic" because some physical quantities, such as opacities or nuclear reaction 
rates, have been tuned to better reproduce the sound-speed profile in the radiative 
zone deduced from the analysis of p-mode observations. One of its aims was to 
improve the prediction of the neutrino flux. This model includes turbulent diffusion 
in the tachocline and the chemical composition of Grevesse and Noels (1993). 

As we want to compare the impact of different physical processes in our work, five 
other solar models have been used. They mainly differ in three physical inputs: the 
microscopic diffusion (including the gravitational settling), the horizontal turbulent 
diffusion in the tachocline, and the abundances of chemical elements. The model 
called no diffusion does not treat any diffusion processes, whereas all of the following 
models include it through the prescription of Michaud and Proffitt (1993). Although 
today this model is obsolete, it is a good test for the detection methodology, since the 
g-mode frequencies are very different as compared to the seismic ones and therefore, 
the analysis should be able to clearly separate them. 

The second process, the tachocline mixing, is introduced in the code following the 
prescription of Spiegel and Zahn (1992). The models including this process are called 
tacho while the ones without this turbulent effect are named std. This turbulent hori- 
zontal motion contributes to the thinness of the tachocline responsible for the mixing 
of chemical elements. Finally, the chemical composition is either taken from Grevesse 
and Noels (1993) (GN93 hereafter) or from the recent release of Asplund, Grevesse, 
and Sauval, 2005 (As05 hereafter). The difference between these two chemical abun- 
dances is a decrease of the heavy element by « 30%. Such differences substantially 
affect, for instance, the Rossland opacities and lead to noticeably different thermal 
stratifications. The models using the old abundances have the suffix 93 whereas the 
models including the new abundances have the suffix 05. Table 1 summarizes the 
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main characteristics of all of the models used. Finally, to use a model calculated 
with another stellar evolution code, the following study is also done with model S 
from Christensen-Dalsgaard et al. (1996), which is a standard model including the 
microscopic diffusion and using GN93. 

Table 1. The main characteristics of the solar evolution models used to 
compute the g-mode frequencies. 



Microscopic Diffusion in 
Solar model diffusion tachocline GN93 As05 



no diffusion x 

std93 x x 

std05 x x 

tacho93 x x x 

tacho05 x x x 



From each stellar structure, eigenfrequencies have been computed with the adi- 
abatic pulsation package (adipack provided by J. Christensen-Dalsgaard, available 
on http://www.phys.au.dk/~jcd/adipack.ri/). The computation of g-mode frequencies 
can become very sensitive to the numerical approach, especially to the treatment 
of the boundary condition at the center. Thereby, the numerical uncertainties on 
the computed g-mode periods become large at very low frequency and we should 
be extremely cautious. However, in the range used for this study (between 2 and 
8.5 hours, see Section 3), it has been shown that, for a given model, frequencies are 
modified by less than 50 nHz by changing the numerics (see Mathur, 2007), leading 
to a maximum uncertainty around 50 seconds on the mode periods. Such errors are 
smaller than the observed deviations between the different models used here. In fact, 
they are too small (by more than an order of magnitude) to be correctly distinguished 
by the methodology developed in this paper. 

It is usually considered that the high-frequency g modes (above 150 /iHz) are 
very sensitive to the physical processes included in the models (the structure) giving 
up to 5 /iHz of difference between the models, while at lower frequencies - inside 
the asymptotic region - the frequency differences are smaller than a fraction of a 
/iHz and they could be more useful to study the dynamics of the solar interior (for 
example see Mathur et at, 2007). Indeed, we have already seen that the difference 
of the periods of two consecutive g modes with the same degree (£) and successive 
radial order (n) is almost constant. In Figure 1 the difference (APi) of the central 
frequencies (m=0 components) of the g modes have been plotted for all the models 
described in Table 1. The asymptotic regime is reached for periods above « five 
hours and the differences between all of the models considered here are within one 
minute (from 24 to 25 minutes). Thus, it is going to be very difficult to distinguish 
between all these models. 

Moreover, until now we have only considered the effect of the structure on the 
g-mode frequencies but we cannot neglect the effect of the rotation of the solar 
interior on the other m-components of the modes. Indeed, in Figure 1 we have also 
drawn the mil components for one of the solar models, the std05, considering a 
core rotating as a rigid body, like the rest of the radiative region. The dispersion in 
Afi introduced by the dynamics (without even considering a rapid rotation in the 
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Figure 1. Separations in period (APi) of the m=0 components of the g modes computed using 
the models described in Table 1. For clarity, we have plotted the mil components only for one 
model: the std05 (with a solid rotation of Q rz = 433 nHz in the radiative zone and a rotation 
in the convection zone that reproduces the inferred profile deduced from helioscismology data) . 



core, just a rigid rotation profile) is already of the same order of the structural effects 
due to the different physical processes and abundances considered. That means that 
both effects cannot be easily disentangled and we have to study both at the same 
time. 

At low frequency, it seems more appropriate to study the g modes in period, 
and not in frequency. Depending on the physical processes taken into account in the 
solar models, there can be important differences between the periods of the g modes. 
Thus, we could expect to use these low-frequency g modes to diagnose some physical 
processes or abundances taken into account in the models. 

To complete this study we need to model the dynamics -the rotation rate- of 
the solar interior, in particular, in the regions where the g modes are trapped, i.e. 
the radiative region. Today, the solar rotation rate is known well down to the solar 
core (Thompson et at, 2003) but it is still uncertain in the deepest layers. We have 
thus computed the g-mode splittings corresponding to ten artificial rotation profiles 
f2(r, 9). These artificial profiles have a differential rotation in the convection zone 
(that mimics the real behavior), and a rigid rotation from 0.7 down to a given 
fractional radius r\ = t/Rq (0.1, 0.15, and 0.2) which is equal to f2 ra( j=433 nHz. 
From each of this internal fractional radius (ri) down to the center we have a step-like 
profile having a rate of three, five, and ten times larger than the rest of the radiative 
zone. Finally, we have also computed a rigid rotation rate down to the center of the 
Sun which is considered as our reference rotation profile. 



3. Methodology 

To compare the signature of the p-mode asymptotic properties between real and 
modeled spectra or even between two different synthetic spectra, we have used the 
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method developed in Garcia et al. (2007). Indeed, a numerical simulation of the 
g-mode power spectrum has been used (see Section S.2 of the on-line material of 
Garcia et al. (2007) for a full description of the simulation). Briefly, the central 
frequencies of the g modes have been computed as described in Section 2. Their 
amplitudes have been considered all the same following the calculations done by 
Kumar, Quataert, and Bahcall (1996) where it has been shown that all the modes 
in the studied regions have roughly the same amplitudes. We have also simulated 
the g-mode amplitudes following different laws (e.g. a linear variation, a parabolic 
decreasing with frequency...). In all the cases, the results were qualitatively the same. 
Finally one-bin line widths have been used but we have checked that the results do 
not significantly change when wider modes are simulated. Since we are only interested 
in the power spectrum, the mode-phase information is lost. It is also important to 
note that all the computed simulations do not contain any noise because otherwise 
we could not quantify between the effects coming from the models and those coming 
from the different noise realizations. Therefore, this study represents what we could 
obtain in the best conditions. The addition of noise (both instrumental and solar) 
would affect the data and could probably degrade the results. 

The starting point of this analysis is the oscillation power spectrum density. We 
will call it hereafter the first spectrum (PS1). The real first spectrum is computed 
from GOLF time series calibrated into velocity (Garcia et al, 2005) with a classical 
periodogram; the synthetic ones are directly built from the g-mode frequencies and 
rotational splittings derived from the considered models (see top panel of Figure 2). 

The PSl is expressed as a function of the period, since the g modes present a 
comb-like structure in period and not in frequency as it is the case for p modes. 
We have then performed a spectral analysis of the first spectrum, limited to the 
range 2-11 hours (i.e. 25-140 /^Hz), to get the second spectrum (hereafter called 
PS2). It is important to note that the considered PS2 is a Complex spectrum, to 
keep information on both amplitude and phase. For practical purposes, as PSl is 
not regularly sampled in period, we cannot use a classical Fast Fourier Transform 
(FFT), we use instead a general method based on sine- wave fittings. 

The asymptotic signature of the I = 1 g modes is detected in |PS2| 2 as a high 
and broad peak around 22 to 25 minutes. Depending on the physics contained on 
the solar model and on the rotation rate in the radiative region, the shape and the 
position of the APi pattern is modified. Thus in the middle panel of Figure 2 we 
show the resultant PS2 of the same seismic model but using two different rotation 
rates in the core, a rigid rotation in the whole radiative zone and a rotation rate five 
times faster than the rest of the radiative zone starting at ri = 0.15. The model 
of the higher rotation rate has more power at 24.5 minutes and widens the APi 
structure. 

Once this signature (a pattern of peaks) has been detected in the PS2, we can 
locate in PSl what contributes to the emergence of the structure in the PS2. To 
do so, we select in PS2 the vicinity around the pattern at APi and around its first 
harmonic at ^APi. We reconstruct the sum of the sine waves corresponding to the 
amplitudes, frequencies and phases measured in PS2 in these limited regions. In other 
words, we perform an inverse Fourier transform of the PS2 after having applied a 
narrow-band filter around what we consider to be the I = 1 <?-mode signature and 
its first harmonic. 

The maxima of these reconstructed waves indicate the location in PSl corre- 
sponding to the positions where we found peaks that have a repetitive pattern. If this 
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Figure 2. Top panel: Power spectral density, PS1, of t = 1 and 2 modes computed using 
the seismic model and two rotation laws in the core: a rigid core (shifted upwards in green) 
and a core rotating five times faster than the rest of the radiative zone at ri=0.15 (continuous 
black lines in all the panels). Middle panel: Power spectral density, PS2, computed from PS1 
between 2 and 11 hours. The shaded regions are those used to reconstruct the signal in PS1. 
Bottom panel: Reconstructed signals in PS1. In the two lower panels, the green dashed lines 
correspond to the case of a rigid core rotation law. 



signature is indeed the one of the dipolar g modes, the maxima of these reconstructed 
waves correspond to the m-components - depending on their splittings - of the 1=1 
g modes. The bottom panel of Figure 2 shows the reconstructed signals of the seismic 
models and the same two rotation rates in the core as in the previous example. At 
low periods, below five hours, the reconstructed waves are nearly the same in both 
cases, but at higher periods, the model with a higher rotation rate in the core follows 
the separation of the two m-components due to the higher splittings and the position 
of the maxima are different compared to the model using the rigid rotation profile. 

Finally we need to quantify the consistency between our models and the GOLF 
observations. Comparing APi is not very informative because, as mentioned by 
Mathur et al. (2007) all of the theoretical computations of APi are in agreement to 
within a minute and, unfortunately, there is not enough resolution to reach such an 
accuracy in PS2. However, a direct comparison of the periods of the modes shows 
clear discrepancies between the different models (see bottom panel of Figure 2). 
Thus, to compare the positions of the modes in PS1, we compute the Pearsons 
correlation coefficients (e.g. see Press et al, 1992) between the reconstructed waves 
of the models and the reference one (Section 4) or with the GOLF data (Section 
5). If the correlation coefficient is clearly positive, which means that the maxima of 
both reconstructed waves coincide generally, that is to say that in both considered 
first spectra (real and synthetic or both synthetic), the g-mode components are 
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overall at the same position. If the modes in one spectrum are shifted compared to 
the modes in the other one, the correlation coefficient decreases to zero and even 
becomes negative. A large negative correlation indicates that each mode component 
in one spectrum is in between two mode components of the other spectrum. 

The range used for the correlation analysis has been restricted between 2 and 
8.5 hours (33 to 140 /zHz). The reason to reduce the low-frequency limit is twofold: 
on one hand, the dynamical effects on the periods are so important when high core 
rotation rates are taken into account that the m-components of different modes 
£ = 1 and consecutive orders (n) are mixed and the interpretation of the results is 
more complicated; on the other hand, the calculation of the frequencies at such rapid 
periods (above nine hours) could be less reliable and numerical effects can play an 
important role as it has been said in the previous section. 



4. Comparison with the Seismic Model 

We have seen in the previous sections that the period of the g modes in the low- 
frequency region (below 150 /iHz) is dominated by the physics and the dynamics of 
the radiative region, and mainly inside the solar core. We start studying the influence 
of the dynamics on the correlation coefficients. To do so, we apply the methodology 
described in the previous section to the seismic model with a rigid rotation profile 
that is taken as our reference model. Then, we compute the correlation coefficient 
between the waves reconstructed for this reference model and the ones for the seismic 
model with different rotation profiles. The results are shown in Figure 3. 

We can note that the higher the rotation rate in the core, the lower the value of 
the correlation will be (i.e. when the difference between the reference and the studied 
model is bigger, the correlation is smaller). For a given rotation rate in the core, if 
we shift the value of r, from 0.1 to 0.2 R@, the correlation decreases a little showing 
that there is low sensitivity to the fractional radius r,;. Then if we had only one 
reliable model with defined physical processes, we could infer some information on 
the rotation profile thanks to this correlation. We would be able to say, for instance, 
that a very rapid rotation in the core is very unlikely to be compatible with our 
reference model (seismic with rigid rotation profile). 

However, another parameter has to be taken into account: the solar structure. In 
Figure 4, we have plotted the correlation coefficients between the reconstructed waves 
for the seismic model with a given rotation profile and another model with the same 
rotation profile. As each point on the plot is calculated with two different models 
but with the same rotation profile, it indicates the sensitivity of the correlation to 
the different physical parameters included in both compared models. 

We can see in Figure 4 that there are two groups of correlations clearly separated. 
The first one contains the models with the old abundances of GN93 that give corre- 
lations above 70% whatever rotation profile to which they are attributed. Whereas 
the second group of models containing the model without diffusion and the models 
that include the new abundances of As05, gives correlations below 30% and even 
negative. Thus the g-mode frequencies extracted from these two groups of modes are 
not compatible among themselves. While this result is already known (Zaatri et al., 
2007; Mathur et al., 2007) it proves the ability of this methodology to separate both 
types of models. Obviously, there are tiny differences between other models such as 
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Figure 3. Correlation coefficients between the reconstructed waves for the seismic model with 
a rigid rotation profile and the ones for the seismic model with other rotation profiles. All of 
the rotation profiles have a differential rotation in the convective zone and a rigid rotation 
from 0.7 down to n (= t/Rq): 0.1 (dotted line), 0.15 (dashed line), and 0.2 (dottcd-dashed 
line). The rotation rate below is constant: 3C ra( j (crosses), 5f2 ra( j (asterisks), and 10f2 ra[ j 
(diamonds), where f2 ra d is the rotation rate in the rest of the radiative zone (433 nHz). 
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Figure 4. Correlation coefficients between the reference model (seismic) with a given rotation 
profile depending on the rate in the region below n as described in Figure 3 (plus the rigid 
rotation profile (triangles)), and another model with the same rotation profile. These models 
(described in Section 2) are: std93 (blue), std05 (red), tacho93 (light green), tacho05 (dark 
green), no diffusion (black), and model S (pink). 
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the std93 and the tacho93 models because the frequency differences between the g 
modes computed from them are very small. 

Moreover, we can see another systematics in Figure 4 as the correlation coefficient 
decreases with higher rotation rates in the core. When higher rotation rates are 
considered, the reconstructed waves have to match both mil components of the 
1=1 modes. Therefore, the structure of the waves is more complicated allowing 
a higher precision in the calculation of the correlation and more sensitivity to the 
positions of the modes (the structure parameters of the model) . 



5. Comparison with the Real Data 

To compare our models with the real data, we use 3481 day-long GOLF time series 
that were also used in Garcia et al. (2007) (see that paper for more details on the 
data analysis process). Figure 5 shows the results of the correlation between the 
reconstructed waves obtained from the GOLF data and those from the models. In 
this figure, instead of having one free parameter, (physical processes or rotation), 
the results depend on both structure and dynamics. However, by mixing them, it is 
possible to extract some information on the data. 
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Figure 5. Correlation coefficients between the reconstructed waves from GOLF data and 
the ones from all of the models with the various rotation profiles. These models (described 
in Section 2) are: std93 (blue), std05 (red), tacho93 (light green), tacho05 (dark green), no 
diffusion (black), model S (pink), and seismic (light blue). 



As in the case of the correlation with the reference model, we see a clear separation 
between the same two groups of models: on one side we have all the models using 
GN93 opacities and, on the other, the models using As05 opacities and the model 
without diffusion. This was an expected result. Indeed, we knew that GOLF is well 
correlated with the seismic model for rotation rates in the core not faster than five 
times - on average - the rotation rate in the rest of the radiative region whatever the 
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considered fractional radius (r,). As we have seen in the previous section, the models 
containing the new abundances and the one without the microscopic diffusion are 
not compatible with the models including the old chemical abundances. Therefore, 
they could not be compatible with the GOLF data. Actually, models that are close 
to the seismic model (that is to say, with the GN93 composition) associated with a 
rotation rate in the core Cl c < 5 f2 rac i present the best compatibility with the data. A 
rotation profile having a rate five times the rotation rate in the rest of the radiative 
zone, gives the highest correlation coefficients. 

It is also important to note that the models using As05 and the one without 
the treatment of the microscopic diffusion have, in most cases, negative correlations. 
The g-mode frequencies of these models are shifted as compared to GOLF and the 
maxima of the reconstructed waves are placed in between those of the GOLF data. 

We know that the largest difference of (/-mode frequencies can reach up to 5 /iHz 
between the models with the GN93 abundances and the other group of models. In 
contrast to that, among the group of models presenting the highest correlations in 
Figure 5, the frequencies are shifted by less than 1 fiKz (see Mathur et al, 2007). 
Therefore, we can say that using this technique on GOLF data, we are able to dis- 
tinguish between structural parameters of the solar models (such as the abundances) 
when the frequency differences induced by these modifications are of the order of 1 
/iHz or even larger. If solar models that differ by various physical processes, produce 
g-mode frequencies with differences lower than a /iHz, the present methodology 
would not be able to distinguish between them. To improve our detection capability 
we should enlarge the analyzed region toward lower frequencies where we have seen 
that the g-mode periods are sensitive to the models. 



6. Conclusion 

In the present paper, we have studied the changes in the signal reconstructed from 
the period separation APi, which vary when the physics of the solar models are 
modified, and we have verified whether it was possible to detect these changes in the 
g-mode asymptotic signature with the present uncertainties. 

By comparing our different models with the reference one, we have learned that 
the methodology presented here offers enough sensitivity to distinguish between two 
sets of models: first, diffusive models using the old GN93 abundances; second, models 
using the new As05 abundances and the one not including microscopic diffusion. 

This work shows that the g-mode signature detected in GOLF data clearly favours 
the models including the old abundances from Grevesse and Noels (1993). In this 
sense, the analysis of the g-mode signature recently detected is in agreement with 
previous results coming from the p-mode analysis showing a larger discrepancy in 
the sound-speed profile when we use the new abundances from Asplund, Grevesse, 
and Sauval (2005). 

However, we have to keep in mind that all of these results do not mean that the 
new abundances are incorrect. They only tell us that the classic models including 
them are not compatible with the observations. But it could be possible to obtain 
a better agreement by changing other physical quantities of the models such as the 
opacities or the reaction rates or by including new processes to compensate for the 
effects of the reduction in the metallicity. 
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